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ABSTRACT 

Among all galactic microlensing events, those involving a passage of the observed source 
star over the caustic created by a binary lens are particularly useful in providing information 
about stellar atmospheres, the dynamics of stellar populations in our own and neighbouring 
galaxies, and the statistical properties of stellar and sub-stellar binaries. This paper presents 
a comprehensive guide for modelling and interpreting the lightcurves obtained in events in- 
volving fold-caustic passages. A new general, consistent, and optimal choice of parameters 
provides a deep understanding of the involved features, avoids numerical difficulties and min- 
imizes correlations between model parameters. While the photometric data of a microlensing 
event around a caustic passage itself do not provide constraints on the characteristics of the 
underlying binary lens and do not allow predictions of the behaviour of other regions of the 
lightcurve, vital constraints can be obtained in an efficient way if these are combined with 
a few simple characteristics of data outside the caustic passages. A corresponding algorithm 
containing some improvements over an earlier approach which takes into account multi-site 
observations is presented and discussed in detail together with the arising parameter con- 
straints paying special attention to the role of source and background fluxes. 

Key words: gravitational lensing - methods: data analysis - stars: atmospheres - binaries: 
general - stars: low-mass, brown dwarfs - Galaxy: stellar content. 



1 INTRODUCTION 

The majority of the galactic microlensing events that have been 
detected or are currently detected are compatible with both the ob- 
served source star and the compact massive lens being approxi- 
mated by point-like objects. However, about 2 to 3 per cent of the 
events involve peaks lasting from a few hours to a few days whose 
shape is characteristic for the source crossing a fold-caustic line 
which is created by a binary (or multiple) lens object. These fold- 
caustic passage events provide valuable information about both the 
source star and the lens object not being extractable from 'ordinary' 
microlensing events. 

Microlensing events compatible with a point-like source at 
distance Ds and a point-like lens of mass M at distance Dl in- 
volve only one parameter that is related to physical properties of 
lens or source and carries a dimension, namely the time-scale of 
motion tn = 9e/h, where /i denotes the proper motion of the 
source relative to the lens and 9e is the angular Einstein radius, 
given by 



0e = 



4GM D s - -Dl 



(1) 



The time-scale of motion £e is therefore a convolution of the lens 
mass M, the lens distance Dl, and the proper motion /i, which 
cannot be measured individually, so that the power of ordinary mi- 



crolensing events for the determination of the mass function and 
phase-space distribution of stellar populations in our own or neigh- 
bouring galaxies that have caused these events is severely lim- 
ited (Mao & Paczvnski 1996). In contrast to this, the observed du- 
ration of a fold-caustic passage t~t provides a second time-scale 
which is related to the angular size of the source star An es- 
timate for 6>* can be obtained from the stellar flux and spectral 
type, so that tjr provides a lower limit to the proper motion n by 
yielding its component = 0*/t* perpendicular to the caus- 
tic, and jj, = /i _L /(sin0) itself follows with a determination of 
the crossing angle qb from modelling the full observed lightcurve. 
The measurement of the relative proper motion between lens and 



source for a microlensing 


event toward the Small Mag 


ellanic Cloud 


(SMC) (Afonso et ai]|l998t IXlbrow et al. 1999c b; 


Afonso et alj 


2000; Alcocketal. 1999 


|Rhieet al| 1999: Udalski et al. 1998) in- 



dicated that the lens in located in the SMC rather than in the Galaxy, 
giving support to the hypothesis that the majority of microlensing 
events toward the S MC are due to self-lensing among its stellar 
populations (e.g. lDi Stefanol2000l) . 

In addition to providing a measurement of the proper motion 
/i, lightcurves around caustic passages are particularly sensitive to 
effects caused by the parallactic motion of the Earth around the Sun 
or by the orbital motion o f the binary lens I Hardv & Walke Jl995t 
iGould & Androno\ll 999). Observations of these effects and the de- 
termination of the corresponding parameters (e.g. Dominikl l 19981) 
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together with the determination of the proper motion n break the 
degeneracy between the lens mass M and its distance Dl yielding 
measurements of the individual quantities (up to possible discrete 
ambig uities). One such measurement has been reported bv lAn et all 
12002). 

The characteristic and distinctive shape of the lightcurve 
around caustic passages allows to obtain stronger constraints on the 
parameters of the binary lens than those resulting from binary lens 
events that show only weak deviations compared to single lenses. 
Therefore fold-caustic passage events provide the main source of 
information about the statistical distributions of the properties of 
stellar and sub-stellar binaries such as the total mass, mass ratio, 
semimajor axis, or orbital period. 

The most valuable results arising from caustic-passage mi- 
crolensing events so far have been on the study of stellar at- 
mospheres. The strong differential magnification that arises as 
the source passes over the caustic allows to resolve surfaces 
of source stars at Galactic distances providing a powerful and 
unique technique to study stellar atmospheres of a variety of dif- 
ferent types of stars ther e by prob ing existi ng theoretical m odels 
ISchneider & Weid Il987t TRhie & Be nnett! Il999t iGaudi & Gould 
Il999h . Measurements of limb-darkening coefficients from binary 
lens m icrolensin g events have be en published for three K-giants 
Albrowetal|l999i! , l200Ct iFields et aljl2003l). one G/K-su bgiant 
Albrowet all2001al) and a Solar-like star <Abe et all2003l) in the 
Galactic bulge as well as for an A-dwarf in the SMC ( Afon so et alJ 
2000), where the source sweeps over a cusp for two of these events. 
In addition to a dense photometric coverage during the course of a 
caustic passage, the observation of a temporal sequence of spec- 
tra pro vides a probe of the spatial distribution of chemical ele- 
ments iGaudi & Gould fl 9991) using the observed variation of re- 
lated spectral lines as indicators. For a galactic bulge microlensing 
event in 2000, several low-resolution spectra have been taken with 
the FORS1 spectrograph at the VLT (ESO Paranal, Chile) show- 
in g a significant variatio n of the equivalent width of the Ha-line 
bv lAlbrowetai]<2001bh . while two high-resolution spectra have 
been taken a t the K eck telescope (Mauna Kea, Hawaii, USA) by 
ICastro et alJ <200ll). and both measurements have later been dis- 
cussed bv lAfonso etail <200ll) . In 2002, a dense sequence of high- 
resolution spectra has been obtained for another event with the 
UVES spectrograph at the VLT where significant equivalent-width 
variations in the m ost prominent lines Ha, H/3, an d Ca II have been 
detected JCassan et all2004lBeaulieu et all2004l) . 

In addition to the events already discussed in the litera- 
ture, a dense high-precision photometric coverage of caustic- 
passage regions of several further microlensing events allowing 
the determination of limb-darkening coefficie nts has been ob- 
tained by PLANET 1 <Albrow et alJ Il998l : iDominik et ail l2002t) 
during recent years, and additional data on caustic-passage 
events have also been collec ted by other micro lensing col- 
laborations, namely E ROS 2 jAfonso et all \200A MACHO 3 
jAlcock et alJl2000Jblcl). OG LE 4 iUdalski. Kubiak & Szvmanskil 
Il997t IWozniak et al.J l200lh . MOA 5 ISumi et alJ 120031) . MPS 6 



jRhie et alj[l999l) . and MicroFUN 7 <Yoo et aljEool . With cur- 
rently ~500 microlensing alerts per year issued by OGLE-III 8 
and ~60 alerts per year issued by MOA 9 , about 10-15 caustic- 
passage events per year can be expected for which limb-darkening 
coefficients of the source star can be obtained by a dense high- 
precision photometric coverage, while, depending on allocated re- 
sources, spectroscopic measurements should be feasible for a few 
events per year. 

The aim of this paper is to provide a comprehensive guide 
for modelling and interpreting microlensing lightcurves near fold 
singularities which are produced by extended sources and binary 
lenses, and for finding all suitable models for microlensing events 
involving fold-caustic passages. In addition, the derivation of phys- 
ical properties of the lens and the source resulting from the mod- 
elling of such events is discussed. This paper outlines the underly- 
ing theory and sets the notation for upcoming analyses of observed 
data. A similar paper about cusp singularities including a compar- 
ison with the relations for folds as discussed in this paper is under 
preparation. 

For lightcurves near a fold-caustic passage, a new consistent, 
general, and optimal set of parameters which are directly related to 
observable features of the observed lightcurve is introduced. This 
choice not only leads to a better understanding of the features of the 
lightcurve by allowing an obvious interpretation in terms of corre- 
sponding properties of lens or source, but also minimizes correla- 
tions between the parameters. In addition, singularities in parame- 
ter space likely to cause numerical problems, e.g. when the source 
size tends to zero or the background flux assumes some special val- 
ues (in particular for flux values arising from image-subtraction), 
are avoided. 

Finding all suitable model parameters corresponding to 
lightcurves consistent with the observed data of a complex mi- 
crolensing event is a non-trivial task due to the large num- 
ber of parameters and intrinsic degeneracies and ambiguities 
between the physical properties of the source and the lens 
system <Mao & Pi Stefanol 1 19951: IDominik & Hirshfeld Il996t 
IPi Stefano & Pernall997t|PominilJl999JbHAlbrow et alJl999bl) . 
Since the lightcurve of sources in the vicinity of the fold is de- 
termined by local properties related to the fold singularity, limb- 
darkening measurements do not require the assessment of models 
for the complete lightcurve along with arising ambiguities and de- 
generacies. However, a model of the full lightcurve is needed for 
determining the proper motion fi, the time-scale of motion £e that 
carries information about lens mass M and distance Dl, and the 
mass ratio q and the separation parameter d = S/9e of the binary 
lens (where 8 is the instantaneous angular separation between its 
components). Moreover, a dense coverage of the lightcurve during 
the caustic passage and in particular target-of-opportunity observa- 
tions of spectra yielding a powerful test of the atmosphere of the 
source star require some prior arrangements pushing the need for 
a prediction of caustic passages. With the lightcurve near a caus- 
tic passage being determined by local properties only, a prediction 
is possible with data taken on the rise to a caustic exit, while the 
data taken over previous caustic passages itself do not provide any 
information, whereas the combination with previous data outside 
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caustic passages can yiel d some constraints lAlbrowet al]|l9 99b: 
Ijaroszvnski & Maol200ll) , 

In general, the correlations between model parameters are 
minimized by choosing parameters that correspond to observable 
characteristics of the observed lightcurve. For binary lens caustic- 
passage microlensing events, the caustic passages themselves show 
some characteristic features. Restricting the full parameter space by 
modelling a caustic passage and performing a search in the lower- 
dimensional subspace holding additional parameters necessary to 
describe the full lightcurve therefore constitutes an efficient method 
for modelling such events lAlbrowet ai][T9 99b). In this paper, an 
improved variant of this approach taking into account multi-site 
and/or multi-band observations is described and discussed in de- 
tail, and the relations between the fold-caustic model parameters, 
the properties of the binary lens, and observable characteristics of 
the full lightcurve such as the time-scale of motion or the source 
and background fluxes are investigated. 

In Sect.|2| the basic theory of gravitational lensing is reviewed 
setting the basic notation used in this paper, while Sect.|5|gives a 
detailed discussion of the fold singularity and its properties. Sect.|4] 
investigates the magnification of source stars in the vicinity of fold 
caustics and Sect.[3]discusses the specific implications arising from 
limb darkening. Sect.|6|deals with lightcurves of source stars in a 
region close to a caustic passage within a microlensing event and 
presents an adequate parametrization of such parts of lightcurves to 
be used for modelling observed data where the model parameters 
directly correspond to observable properties of the lightcurve and 
can be easily understood. The determination of a complete set of 
parameters characterizing the full lightcurve by making use of the 
parameters locally characterizing the caustic region forms the con- 
tent of Sect. Q where two versions of the basic algorithm (differ- 
ing in the consideration of the temporal variation of the magnifica- 
tion of non-critical images) are discussed, followed by discussions 
of constraints arising from non-negative background flux and from 
simple characteristics of the data outside the caustic-passage region 
such as the baseline flux or the time-scale of motion. In Sect. [8] 
the predictive power of the data in the caustic passage region and 
its combination with characteristics of data outside it with respect 
to the determination of the characteristics of the binary lens and 
the perspective for predicting future caustic passages is discussed. 
Sect. [5] finally provides a summary of the paper and its results. A 
table of the used symbols has been attached at the end of the paper. 



2 THE LENS MAPPING 

Gravitational lensing is understood as the bending of light emitted 
by a source caused by the gravitational field of intervening mat- 
ter. Due to this effect, a source located at a position y on the sky 
(in arbitrarily scaled coordinates) will be observed at one or more 
image positions x^K Moreover, the luminosity of the observed im- 
ages differs from the intrinsic luminosity of the source object. In 
general, a lens mapping can be described by the Fermat-Potential 



$>{n,y) = ^ ( x - y) 2 - ipi x ) 



(2) 



lSchneidedll985h corresponding to the arrival time of (hypothet- 
ical) light rays at the observer. Fermat's principle determines the 
actual light rays to satisfy 



y = x - a{x) . 
where 



(4) 



a(x) = V x iP{x). (5) 

Let indices to $ denote its n-th partial derivatives with respect to 
a?ij . . . Xi n , ik being the coordinate index of the fc-th derivative, i.e. 

d d 



(6) 



and let t/>ii,...,i„ denote the corresponding derivatives of if). For (lo- 
cally) continuous derivatives, the order of derivation and therefore 
the order of indices to $ or tf> is irrelevant (Schwarz's theorem). 
For the second derivatives of $, one obtains the relation 



$i 3 (as) = S i3 - ipij (x) = , 
where 8ij denotes the Kronecker symbol 



5ij 



for i 7^ j 

1 for i — j 



(V) 



(8) 



Eq. implies that all derivates of y with respect to Xi k can be 
expressed by means of the Fermat potential $ and its derivatives. 
In particular, the components of the Jacobian matrix J read 



9i 



and its determinant becomes 
det J = $11 $22 - ($12) 2 ■ 



(9) 



(10) 



The Jacobian determinant yields the (signed) magnification of 



an image at x^ l > as 



(ID 



detJ(xM) 

and the total magnification A of a source at y is obtained by sum- 
ming over the absolute values of the individual magnifications of 
its N images, i.e. 



N 

E 



(12) 



(3) 



where x^ and y satisfy the lens equation. 

The lens mapping becomes singular at points a;' crlt ' in 
image space for which the Jacobian determinant vanishes 
(det J(x^ crlt - ) ) = 0) which themselves are called critical points. 
The set of critical points forms critical curves which are mapped 
onto the caustics by the lens equation. The number of images asso- 
ciated with a given source position changes (by multiples of two) 
with the source position if and only if the source crosses a caus- 
tic. Therefore, the caustics divide the two-dimensional space of 
source positions into regions with a fixed number of images, where 
a source is defined to be 'inside' a caustic if it belongs to the region 
with the larger number of images, and 'outside' the caustic other- 
wise. The singularities of the lens mapping can be categorized into 
different types showing their own specific characteristics. 



3 THE FOLD SINGULARITY 

The lowest-order singularity i s the fold, for whic h the following 
conditions need to be fulfilled 1 Schneider & W eill 19921) : 



which relates source and image positions by the lens equation 



• exactly one of the eigenvalues of the Jacobian matrix is zero, 
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• the tangent vector of the critical curve is not an eigenvector of 
the Jacobian matrix belonging to eigenvalue zero. 



For critical points, the eigenvalues of J are given by 
A (1) =$11 +$22 =TrJ, A (2) =0. 



(13) 



The first condition for a fold then implies that $n + $22 7^ 0. 
Normalized eigenvectors that correspond to these eigenvalues can 
be written in the form 



(14) 



e« = | 


cos 8 \ 
sine i 


, e (2) =( 


— sin 8 \ 








cos 8 J 



where 



= e($i2 $2 



<i> 2 



$11 + $22 



cos U = 



$11 



$11 + $2 



with 

e(z) 



1 for z > 
— 1 for z < 



(15) 



(16) 



The vanishing of the Jacobian determinant yields $11 $22 = 
($i2) 2 > 0, so that |$n +$22| = |$ll| + |$22|. 

Let Xf denote a point on the critical curve and y l — Xf — 
ct(xf) denote the corresponding source coordinate on the fold 
caustic. The eigenvectors of J given by Eq. 1 141 then constitute 
a right-handed basis of a coordinate system defined by 



x(t) = x f +1Zx'(t), 
y(t) = y f + TZy'(t), 



where 

n = 



cos 8 
sin 8 



— sin 8 
cos 8 



(17) 
(18) 



(19) 



so that 8 is the orientation of the (xi,x' 2 )- and (2/1,3/2)- coordi- 
nates with respect to the (asi, £2)- and (3/1, ^-coordinates, re- 
spectively. Since every multiple of an eigenvector is an eigenvector 
itself, — (e' 1 ', e' 2 ') also constitutes a right-handed basis being re- 
lated to (e (1 \ e^ ) by a 180° -rotation, equivalent to replacing 8 
by 7t + 8. Eq. (T71 implies that 

dxi 
dxT ' 



1Z 



(20) 



The Jacobian matrices J(£i, X2) and J'(x'i, x' 2 ), where Jy = 
dyi/dxj and Jy = dy^/dxj, are related by 

J' = 1l T J1Z, (21) 

so that J' becomes diagonal. In analogy to Eq. ^6), let indices to $' 
denote derivatives with respect to x\. The invariance of trace and 
determinant under the transformation given by Eq. <2U yields 

$'n = A (1) = $11 + $22 , $22 = A (2) = . (22) 

In the chosen coordinates, the gradient of the Jacobian deter- 
minant reads 



(23) 



so that a tangent vector to the critical curve has to be proportional 
to (— $2??> $i??) T . The refore, the conditions for a fold become 
fcchne ider & W eill 199 A 



• $'n ^ 0, $22 = 0, 

• $222 + 0, 

which guarantee that the gradient of the Jacobian determinant does 
not vanish. Taking into account these conditions and the specific 
ch oice of coordinates, the lens equation near a fold can be expanded 
as iSchneideret all 1991 p. 183f.) 



y[ = &ux'i + - $'122(^2) +$'112^1^2, 



(24) 

2/2 = ^^ii2(Xir + t>i 22 x 1 x 2 + -<f 222 (x 2 )- . (25) 
Using this expansion, the Jacobian matrix near the fold reads 

j/ = ( $ii + $'112^2 W + $'1222:2 \ (26) 

V $il22:i + $'122^2 $'l22»'l + $222^2 /' 

and to lowest order its determinant becomes 

det J' = $'u($i 12 xi + $222x2) , (27) 

in accordance with Eq. J23t . The condition det J'(cc'' cllt ') = 



yields the critical curve as $'112 x'[ crit) + $'222 x'^"' = 0, so that 
with Eqs. 1241 and 1251 . the caustic is obtained (to lowest order) as 
the parabola 

2 

(28) 



/(crit) 



,(crit) $'ll2$222 

y 2 = 



($122)^ / /(crit) ,2 



2($il) 2 $222 

For source positions along its symmetry axis (y[ — 0), Eqs. J24I 
and J25j yield two images 



x 2 = 



$'12 



V2 ■ 



2y' 2 



(29) 



(30) 



for y' 2 > which degenerate into a single (critical image) for y' 2 = 

0, while there are no images for y' 2 < 0. It follows that the normal 
to the caustic pointing to the inside is in the positive y' 2 -direction, 

1. e. rif = (0, 1) T , so that in (ja, ^-coordinates, tit = e' 2 ' is the 
inside-pointing caustic normal at y { . Neglecting y 2 against y/y^, 
Eqs. J2Z). I2D, and GHJ yield the magnification due to the critical 
images as 



Acrit(0,^) = W^G(^), 
V Vi 



where 

Ri = 



> 



(31) 



(32) 



Wi) 2 1*222 1 

measures the strength of the strength of the fold caustic and can 
also be interpreted as a characteristic distance scale. 

The relation given by Eq. 1201 implies that derivatives of the 
Fermat-potential $ with respect to x\ and Xi are related by 

2 

' $h TZ kt TZi 3 , (33) 



E 

k,l=l 



$' 



I ,m.n— 1 



(34) 



the first of the equations being equivalent with Eq. \2\\ . Explicitly, 
the most relevant derivatives for a fold caustic read 



$'11 



cos 2 8 $11 + 2 sin 8 cos 8 $12 + sin 2 8 $22 , (35) 
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cos 3 9 $111 + 3 sin 2 9 cos 9 $112 - 
— 3 sin 9 cos 2 9 $122 + cos 3 9 $222 



(36) 



By inserting the explicit values for sin 9 and cos 9 given by 
Eq. {B}, one reveals $' n = $n + $22, Eq. <22l . 



4 MAGNIFICATION OF SOURCE STARS 

The magnification of a point source in the vicinity of a fold caustic 
can be decomposed into the contributions by the two critical images 
and the remaining images, i.e. 



(37) 



If one neglects the curvature of the caustic and any changes of the 
lens properties in the direction perpendicular to it, Eq. 13 1 1 together 
with the relation det J = det J' yields the magnification for a point 
source at y due to the critical images as 

/ \-l/2 

^t(vx;ft)=fgj O(yx), (38) 

where y ±_ = (y — y f ) ■ nt is the distanc e perpendicular to the caus- 
tic (c.f. IScnneider & Wei3ll98d ISchneider et alJll99i p.l86f; 
lAlbrow et alJll99 9b). For the validit y and limits of t his approxi- 
mation, the reader is referred to lGaudi & Pettersl 120021) . 

Let us consider a circular source star with radius p* and sur- 
face brightness /' s '(p) = I C' s '(p) f° r a n l ter s as function of 
the fractional radius ^ p ^ 1, where / denotes the average 
value, and £} s \p) is a dimensionless function describing the ra- 
dial stellar brightness profile, which is normalized to yield unity if 
integrated over the source, i.e. 



(39) 



With i"' s '(p) being composed of the specific intensities I x (p) for 
the wavelength A by means of the transmission function B^ 3 \\) 



as 



/ w (p) = / I (p) B [S '(X) dX , 




(40) 



and I x (p) = I X £ A (p), where / A denotes the average intensity and 
£ A (p) a dimensionless normalized profile function for the specific 
wavelength, one obtains 



j i x B (s \x)e( P )d\ 



e s, (p) 



f J A BW(A)dA 



(41) 



For the source being centered at y, one obtains the magnifi- 
cation by convolving Eq. <38t with the brightness profile £ (p), 
yielding 

A cvit (y ± ,p*;Rr,Z (s) ) = A° CTit (yx/p* , p*/R { ; C W ) 
1/2 



(42) 



where Gf is a universal fold-caustic profile function depend- 
ing on the brightness profile £^ s ' of the source given by 



ISchneider & Weilll987h 



max ( 1 — 7} , 1 ) 



71 



max(l — 77, — 1) 



v 7 ^ 



C (s) (tJpI+pI) dP|| dp, 



(43) 



In contrast to several previous papers, r\ is defined so that the source 
is completely outside the caustic for r\ < 0, where G\ ^ 3 ') = 0, 
completely inside the caustic for 77 > 2, and its center crosses 
the caustic for rj = 1. As can be seen directly from Eq. 1421 . 
^4crit (y± , p* ; Rf ; £ ) depends only on the two ratios y± / p, and 
Pi,/R{. Explicit forms of Gf for selected brightness profiles £' s ' 
will be discussed in Sect. [5] 

Compared to A£ rit , the magnification due to the remaining, 
non-critical images is varying slowly in the vicinity of a caustic, so 
that it can be expanded as 



A p othcl (y)=A t + (y-y f )-(VA) f , 



(44) 



VA)f 7^ 0, while otherwise the use of higher-order 



where At - 
sumed that 

terms in the expansion would be indicated. 10 Both At and (V-A)f 
can be expressed by means of derivatives of the Fermat potential 
«J> evaluated at the N n on-critical images x ^ of the source at y f . 
As outlined in detail bv lWitt & Mad (l995), the lens equation for 
a binary lens can be written as 5th-order complex polynomial in 
z — xi + 1x2 which can be solved numerically for the the images 



x^ + 1X2' by standard root-finding routines. While At is 



simply given by 

JV 

Al = y- 1 

|$ n $22 - ($12) 2 

i — 1 



(45) 



the gradient (Vi){ can be determined as follows. 

For an image at xo and the corresponding source at y Q , deriva- 
tives with respect to the components of x and y are related by 



V„ = J 1 (x )V a! , 
so that (V.A)f reads 

N 

(VA) f = ^J- 1 ( a; (i) )V*K* W 

i=l 
N 

Esign(det J (a; 
jdetJ(a;( l ) 



(46) 



1 = 1 

£ sign(det J(^-))) j_ 1(x(i)) Vx ^ J(x(i)) (47) 



With 



and 



det J 



$22 
_$ 12 



_$ 12 
$11 



(48) 



10 The condition (VA)f 7^ is thought to hold for all fold singularities 
of binary lenses, since (VA)j = would require a symmetry of Af to 
spatial variations in directions both parallel and perpendicular to the line 
connecting the two lens objects, which should only be achievable for source 
positions in between the lens objects, where however a saddlepoint rather 
than a maximum or minimum is expected. 
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V x det J 



$11*122 - 2 $12*112 + $22*111 
$11$222 — 2 $12$122 + $22$112 



(49) 



(V.A)f expressed by means of the derivatives of the Fermat poten- 
tial $ becomes 
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-$11$22$222 + 

+ [$ll$22 + 2($l 2 ) 2 ]$122- 

- 3 $ 12 $22$112 + $22$111 
$11 $222 - 3 $11 $12 $122 + 

+ [$ll$22+2($i 2 ) 2 ]$112- 

- $12$22$111 



(50) 



/ 



Due to symmetry, the magnification of an extended circular 
source A ot h CI (y, p + ; £^ s '), obtained by integration of the convolu- 
tion of the magnification ^4-o thor as given by Eq. J44> with the radial 
source brightness profile, turns out to be equal to that of a point 
source at its center, i.e. A othcl (y, p*;£ (s) ) = Al thel (y) unless the 
source size exceeds the range for which the expansion is a fair ap- 
proximation, so that the total magnification of an extended circular 
source reads A iold (y, p*; £ (s) ) = A cilt (y, p*; £ (s) ) 



•^otherCl/)- 



5 LIMB DARKENING 

With # denoting the angle between the normal to the stellar sur- 
face and the direction to the observer, so that cos i3 — yj\ — p 2 , 
normalized brightness profile functions 

t( s ) („\ — (1 _1_ irt l e )\ fl _ „ 2 W 2 



(p) = (l+p/2)(l-p- 
are proportional to cos p #, and their linear superposition 



£ w (p ; r« } ...r 



(a) 

{P kl } 



= 1 



E^ tew 



(51) 



,(52) 



where pi > 0, provides a popular class of models for limb dark- 
ening. The contribution of individual base profiles £ j (p) to the 
full brightness profile function £ < - s ' (p) is measured by the k\ limb- 

( s) 

darkening coefficients F r , ^ 1, which moreover have to 
fulfill the condition 

fei 



E ] 



{Pi} 



^ 1. 



(53) 



The stellar brightness profiles for uniform brightness 
(p = 0), square-root (p = 1/2), linear (p = 1) and quadratic 
limb darkening (p = 2) as a function of the fractional radius p are 
shown in Fig. Q In general, larger powers of p provide stronger 
limb darkening. 

By inserting the brightness profile mto Eq. 1431 . 

one obtains its fold-caustic magnification function Gf ^ p y (rj) = 
Gf{m t\ a J } ) as lSchneider&Weilll987l) 



G{,{ p }(v) = —j= 



1 (1 + !)! 



max(l — 77,1) 



(i 



^Jx + rj - 1 



dx . (54) 



max(l — 77, — 1) 

For even p (including p — 0), the fold-caustic magnification 
function Gf r p \ can be expressed by means of the complete ellipti- 
cal integrals of the first and second kind, K(x) and E(x), respec- 
tively (e.g. iGradshtevn & Rvzhikll994T) . while Gf r p j becomes an 
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Figure 1. Stellar brightness profiles , (p) 
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(l+p/2)(l 

a function of the fractional radius p for source stars with uniform brightness 
(p = 0), square-root (p = 1/2), linear (p = 1), and quadratic (p = 2) limb 
darkening. 



analytical fu nction for odd p. For unif ormly bright sources (p = 0), 
one obtains iScr,neirter& Weill! 9871) 

'0 for rj sC 

for < rj < 2 , (55) 



-(i-n)E 



for rj > 2 



while for linear li mb darke ning (p = 1) the evaluation of Eq. <54t 
yields iSchneider & WagoneJl987l) 



Gt{i}(v) = 



,,3/2 



5 (5 - 2t?) r," 
| [(5-2r ? )r ? 3 / 2 



for 



for rj < 

< rj <: 2 



+ 



(56) 



+ (1 + 27?) fa - 2) 3/2 ] for 77 > 2 



The fold-caustic magnification functions Gf for uniformly 
bright sources (p = 0), square-root (p = 1/2), linear (p = 1) and 
quadratic limb darkening (p = 2) are shown in Fig.|2| With the stel- 
lar surface brightness vanishing at the limb for all power-law pro- 
files with p > 0, the slope of Gf ^ p y becomes zero when the stellar 
limb touches the caustic from the outside (rj — 0), In contrast, there 
is a slope discontinuity at this phase for p = and therefore for all 
brightness profiles that involve a non-zero constant brightness term. 
As the source enters the caustic, the magnification rises to a peak 
and thereafter falls asymptotically with the inverse square root of 
its perpendicular distance to the caustic. The position of the peak 
depends on the brightness profile and occurs for the source center 
being inside the caustic. With more light being concentrated near 
the source center, stronger limb darkening can account for fold- 
caustic peak magnifications similar to that for less limb-darkened 
but smaller sources. On entering the caustic, stronger limb dark- 
ening produces less steep initial rises of the magnification, while 
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Figure 2. Universal fold-caustic magnification function Gj , , (rf) for 
selected source brightness profiles following power-laws in cos $ = 
y/l - p 2 of the form f ^ (p) = (f + p/2) (f - p 2 )^/ 2 : Uniform bright- 
ness (p = 0), square-root (p = 1/2), linear (p = 1), and quadratic (p = 2) 
limb darkening. 



steeper rises occur before a narrower peak at higher magnification 
is reached for a smaller distance of the source center from the caus- 
tic. After the peak, the magnification for extended sources exceeds 
that of a point source, where the difference increases for weaker 
limb darkening. 

For more details about the fo ld-caustic magnification func- 
tion Gf { p } , the reader is referred to lRhie & Bennett! Il999l) . while 
more about the extraction of limb-darkening parameters from the 
observed data during a fold-caustic passage event can be found in 
two recent papers bv lDominTk! i200 4albl) . 



6 LIGHTCURVES NEAR FOLDS 

Let the center of the source cross the caustic at time if and let t| de- 
note the point of time when the source begins to enter or finishes to 
exit the caustic. If the source brightness does not vanish at the stel- 
lar limb, Gf{rj\ S;^) involves a slope discontinuity at r\ = 0, pro- 
ducing a characteristic feature in the lightcurve at time if , whereas 
if is much less easily recognizable from the lightcurve. The choice 
of if as reference point rather than if constitut es an import ance 
difference to the discussion presented bv lAlbrow etafl f 1 999b ) and 
is the key point for avoiding strong correlations between model pa- 
rameters. With t j~ = ± (if — if ) ^ (throughout the paper, the up- 
per sign will refer to a caustic entry, while the lower sign will refer 
to a caustic exit), the source crosses the caustic during the timespan 
2 i* . Let < (j> < 71 denote the angle from the caustic tangent to 
the source trajectory, and iE > denote the timespan in which the 
source moves by a unit distance. During the time ig = iE/(sin <f>), 
the source therefore moves by a unit distance perpendicular to the 
caustic, and tjr = p+t^, so that 



i-if i-if (j_t-tt , 



(57) 



With t denoting an arbitrarily chosen unit time, and 



tt = Rf is > 



(58) 



a characteristic caustic rise time, one can define a caustic rise pa- 
rameter 



> o. 



(59) 



With these definitions, Eq. 1421 yields the magnification due to the 
critical images as 



A crit (t) 



A CI it ( i ■ 



= Cf 



ti 

1/2 



1 **,-£« 



f,7 [ ± *_*!.;£« 



(60) 



The singularity in this expression for t* = (point source) 
can be avoided by defining a function 



Gf (y±,p*;£, 



00) 



(*r 1/a x 

xcr (yl/p*;Z (s) ) 

for pi 

(yl)' 1/2 6(01) 

for pi = 



(61) 



which is continuous in t* = 0, i.e. 



Km <5?I±**L 
tj— o V t t 



= Gf 



t-tt 



t 



-,0;£ 



{■>) 



(62) 



and allows the magnification due to the critical images to be ex- 
pressed as 



Arit(i) = Cf Gf 



t-tr t 



t 



(63) 



(s) 

Consider n lightcurves being observed and let F§ > and 
F^ denote the fluxes of the source and the background, respec- 
tively, for the s-th lightcurve, so that the blend ratio is given by 



9 



00 



Fg /Fg ' and the baseline flux is F l 



F (s) + F 



In 



analogy to the spatial expansion of the magnification v4o thcr due to 
the non-critical images as given by Eq. J44> . its temporal expansion 
around if reads 

Al thel (t)^Al + (t-t* t )Al , (64) 

where A* { = ^o thcr (if*) and Af = ^ ther (if ). The relations be- 
tween the coefficients Af, Af, At, and (VA)f will be discussed In 
SectionQ 

Let us define 



Wf 



F 



,(«) 



x at 



00 



(65) 
(66) 

(67) 



Using the above definitions, the total flux = 



F^ a) [Acrit(i) + ^other + F b' for the s " th lightcurve takes the 



form 
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±&t(t-it) 



00 



(68) 



A fit to the n observed lightcurves therefore involves the 3 + 2n 
parameters t$,ti > 0, F r {s) > 0, F* (s) , and uf. The fluxes F r {s) 
measure the asymptotic behaviour of the lightcurve for source po- 
sitions towards the inside of the caustic (i.e. for t —> ±oo) and can 
be interpreted as the flux of the two critical images that are pro- 
duced at the time t = tf ± i = tf ± (t * + i ) if the extended source 
is replaced by a point source at its center, where the chosen point 
of time corresponds to a unit time after the source center enters or 
before it exits the caustic which itself takes place at t = ff ± t* . 
The parameter describes the temporal variation of A^ ther , where 
± uif measures the rate of the corresponding change of flux in units 
of F r with positive time. Since the change of caustic properties 
in the direction parallel to it has been neglected, the lightcurve is 
affected by the transverse motion only. The observable width of the 
caustic passage given by t* is the product of the source size param- 
eter p, and the time-scale of transverse motion t E , whereas neither 
of the two latter quantities are observables themselves. 

The meaning of the model parameters tj,t+, F^ a \ F f * , and 
Wf is illustrated in Fig. [3] For observed lightcurves, initial guesses 
for these parameters can easily be obtained from recognizable fea- 
tures. The point of time tf is indicated by the arising slope dis- 
continuity, and _F*' S ' is obtained as the flux at tf . The duration of 
the caustic passage 2 t^r extends approximately between the slope 
discontinuity at tf and the change of sign of the curvature on the 
other side of the caustic peak. With two fluxes F| 



and F 2 (s) = F fold 



(£2) taken at two points of time fi and £2 in the 
region where the source is located inside the caustic and can be 
fairly approximated by a point source and the rate of change of flux 
F f of the non-critical images at t { (obtained as the limit of the 
tangent for t—*-tf from the outside), one obtains 



F' 



00 




(F 3 M - J? w )(*!-*?) 



(»)\ 



(69) 



Finally, the slope of any chosen lightcurve k for the source just 



outside the caustic denoted by F f yields ui { 



± F *(k) /F (k)_ 



In the absence of blending, Fg = f° r al l s > so mat Eqs. J65I 



> | M ') (s) 

and 1671 yield a constant ratio between F { and F r , namely 



9t 



F± 



Al 

Cf 



> 0. 



(70) 



In this case, modelling the lightcurve to the observed data involves 
4 + n independent parameters, which can be chosen as tf , tj~ > 0, 
g* { > 0, F* (s) > 0, andw?. 

While in the case of difference-imaging, the flux F B s) repre- 
sents the total background and can have either sign, it is restricted 



to F B S ' ^ for standard photometry. The conditions Fg S; > and 
> 1 imply F f ^ d (t) > for all t in this case, so that Eq. <67l 
requires F f *' s ' > 0. Therefore, a model parameter _F*' fc ) ^ for 
the fc-th lightcurve implies a negative background flux F B . 
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Figure 3. Typical lightcurve for a uniformly bright source entering a fold 
caustic (bold solid line) which illustrates the physical meaning of the model 



parameters t 



> 0, F r (s) > 0, F f 



:(,) 



and uii . For this example, 



with Cj 



2/3 1 and a) f * = 0.08 have been chosen. Also shown are lightcurves 
: for the uniformly bright extended source (light solid line) and 



the corresponding point source (dotted line). The long-dashed line shows 
the contribution by the variation of A^ h (t). 



For sufficiently small ± u>f (f — tf ), F^ d as given by Eq. J68) 
violates the condition F^ d (t) > 0. Although the approxima- 
tion of Mother (t) by its expansion around t — tf only holds for 
|t — tf I <C I Af / Af I , evaluations outside this region may be at- 
tempted during a parameter search. Negative fluxes prohibit the 
calculation of a corresponding magnitude which can however be 
avoided by using a modified expansion of involving an ex- 

ponential function rather than a linear function, which is a purely 
technical modification and has negligible effect on the lightcurve in 
the vicinity of the caustic passage. Let us therefore expand A£ tlier 
as 



A P oth Jt) * ^ + at ( exp i ^lSLJLl > _ , 



Off 



(71) 



where at > 0, so that the local properties A^ thcr (tf ) = Af and 
^hor(*f) =M are fulfilled. 

In analogy to Eq. 1701 . one can define 



f*00 A* , ( 



F, 



Cf 



and 



fff.min = mm 



>0, 



(72) 



(73) 



which corresponds to the lightcurve involving the smallest blend 
ratio. 

With at — CtSimm ar, d &t as defined by Eq. J66I . the total 
flux for the s-th lightcurve reads 



M<0 



G f * ± 



t - tf ti 



{■>) 
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+ flf.min exp i ± 



Wf (t - tf ) 



fff.min 

With the corresponding magnitude being defined as 

(s) f ,\ *(s) „ r , -^fold(^) 



•(«) 



F, 



'(.*) 



(74) 



(75) 



where, in analogy to F f , ?7if denotes the magnitude at time 
tf , the magnitude near a fold-caustic passage takes the form 



(i) = mj w - 
-2.5 lgjl+- 



' ~ i - ' T 



+ 



,min 



(76) 



and lightcurve models involve the 3 + 2n parameters if , > 0, 
<7 f > 0, m f and d) f which reduce to the 4 + n parameters t| , 
i* > 0, (ft > 0, m f and lu { in the absence of blending. 



7 FULL LIGHTCURVE AND BINARY LENS MODEL 
7.1 Constraining the parameter space 

For a galactic microlensing event involving an extended source at 
distance Ds and a binary lens with total mass M at distance Dl, 
let us choose the angular Einstein radius 8e, defined by Eq. Q. as 
the unitlength of our coordinates (xi , X2) and (j/i ,3/2). Usually, the 
resulting lightcurve is characterized by the 7+2n model parameters 
(d, q, Uo, a, to, £e, p*, Fg , F^ ) being described in the following. 
The binary lens itself is characterized by the angular separation of 
its components 8 = d 9e and their mass ratio q. The location of 
the trajectory of the source relative to the lens is described by the 
closest angular separation 80 = «o &e between source and center 
of mass of the lens system and the orientation angle a relative to 
the line connecting the two lens components. The closest approach 
takes place at time to and within the timespan t& = 8e/p the 
source moves by 8e relative to the lens on the sky. The angular 
radius of the source is = p* 8e, and as defined in Sect.|5| Fg S ' 
and F^ denote the sou rce and ba c kgroun d fluxes, respectively. 

As pointed out bv iDominikl ll999bh . the parameters d and 
q are usually strongly correlated resulting in a partial parameter 
degeneracy which can be avoided by instead using the parame- 
ter pair (7, g) for wide binary lenses, where 7 denotes the shear 
and the mass ratio q measures the deviation from a Chang-Refsdal 
lens, or the parameter pair (Q, q) for close binary lenses where 
Q denotes the absolute value of either eigenvalue of the traceless 
quadrupole moment and the mass ratio q measures the deviation 
from a quadrupole lens. Any (d, q) in the following text may there- 
fore also be understood as (7, q) or (Q, q). 

Rather than by the full set of 7 + 2n model parameters 

( s) ( s) 

(d, q, uo, a, to, £e, p*, F§ , F^ ), the lightcurve in the 
vicinity of a fold caustic is completely characterized by 3 + 
2n parameters (if ,t~^, Ft , Ff " , &f ) or 2 + 2n parameters 
(if , t^~, Fx"\ Ff ' s ' 1 ) if one disregards the parameter ujf describ- 
ing the temporal variation of the magnification due to non-critical 
images 11 , as shown in the previous section. Already a comparison 



of the number of parameters shows that the data over the caus- 
tic passage cannot provide an unique model for the full lightcurve 
of the event. Instead, they turn out to be insensitive to the caustic 
crossing angle <f> and the caustic rise parameter £f, as well as to the 
caustic strength Rf and the magnification Af at the beginning of 
the caustic entry or the end of the caustic exit, the last two param- 
eters given by the choice of the singularity y f for a specific binary 
lens characterized by (d, q). Therefore, modelling the caustic pas- 
sage data does not provide measurements of the time-scale tE, the 
source and background fluxes Fg and Fg , the source size pa- 
rameter p±, and of the characteristics of the binary lens, namely the 
separation parameter d and the mass ratio q. 

As will be discussed later in this section, there are however 
some relations between and restrictions on these parameters. In par- 
ticular, the different rise fluxes _F r provide relations between the 

( s) 

source fluxes F^ , and together with the caustic passage fluxes 
Ff yield relations between the background fluxes F^\ while 
the parameter iSf constrains the caustic crossing angle <j> and relates 
it to the caustic rise parameter £f. Moreover, additional constraints 
apply if measurements of properties of the full lightcurve such as 
the time-scale tE, source, background or baseline fluxes F^ , F^ , 
or Full , or blend ratios q^ are taken into account as discussed in 

base' 3 

Sect.1731 

Nevertheless, with a determination of the angular radius of the 
source 8*, e.g. from spectral typing, the caustic passage data yield 
the proper motion perpendicular to the caustic with i* as 



/'• 



ti 



(77) 



whereas the measurement of the full proper motion fj, requires the 
determination of the caustic crossing angle (j> 



/'• : 



t* sin <^> 



/' 



(78) 



where the caustic passage data only yield a lower limit. 

As pointed out by Alb row et alJ Jl999bl) . rather than perform- 
ing a parameter search on the 7 + 2n-dimensional full parameter set 
(d, q, uo, a, to, tE, p+, Fg a \ -Fg 3 '), all suitable binary lens mod- 
els can efficiently be found by using the constraints arising from 
modelling the data around the caustic passage yielding the 3 + 2n 
parameters (if , tjr , F}^ , Ff ^ ,uif), so that only a 4-dimensional 
(or 5-dimensional if u)f is ignored) remaining parameter subspace 
needs to be searched for solutions. This search can be performed by 
creating a suitable parameter grid, and by assessing the goodness- 
of-fit of the arising trial models to the data. For promising regions 
of the parameter grid, the grid may be refined in order to succeed 
towards optimal models and/or the trial models can be used as ini- 
tial guesses for fits in the full parameter space. 

A convenient parametrization of the remaining subspace is 
given by (d, q, £, cf>, ft), where the separation parameter d and the 
mass ratio q are the binary lens characteristics, t stands for the path 
length along the caustic and yields the location of the fold singu- 
larity y t , while the caustic crossing angle <j> and the caustic rise 
parameter Q characterize the source trajectory. 

For a given fold singularity at y f , the local characteristics R{ , 
rif , Af and ( V^4)f are functions of derivatives of the Fermat poten- 
tial of the lens mapping at the critical image Xf or the non-critical 



11 While the inclusion of u)l! as a model parameter yields more adequate 



estimates for other model parameters, its value itself is rather uncertain and 
unreliable. 
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images x^' which are given by Eqs. 1321 . 1 141 . 1451 . and <50> . re- 
spectively, where nf = e' 2 '. In principle, the binary lens model 
yields j4 ot her(t) for any source trajectory and source size. How- 
ever, as pointed out in Sect. |4] within the validity of the approx- 
imation given by Eq. J64> . A ot ^ CI (t) = A ^ thcr (t) which itself is 
simply characterized by the coefficients Af and A] . If the value of 
u>{ is disregarded corresponding to a constant A^ thcl , one obtains 
Af = A{ , whereas otherwise Af = Af — Cf uif t+ . 

7.2 Parameter search disregarding w| 

Let us now investigate how the choice of (d,q,£,(f>,Q) together 
with the 2+2n fold-caustic model parameters (tf , t* , Fr , F f * ) 

Is) (s) 

yields the 7 + 2n parameters (d, q, uo, a, to, is, P*, Fg , Fg ) 
which characterize the full lightcurve. 

Regardless of the caustic crossing angle <f>, the choice of Cf 
dictates the values of all source fluxes Fg S ' and for given (d, q,£) 
the values of the time-scale of perpendicular motion f B , the source 
size parameter = #*/#e, and all background fluxes F B . 

Explicitly, £g is determined by Cf and Ff as 



c£j 
Ff : 



(79) 



while the source size parameter follows from these parameters to- 
gether with t* as 

(80) 



P * t± ^2 - 



'-E 



Together with the conditions Fg > and Af > 1, Eq. 1671 
yields the relations 



*00 _ fO) 
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00 
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(81) 
(82) 



From these equations and the definition of F r , Eq. 1651 . the 
source fluxes Fg S ' and the background fluxes F B S ' of all n 
lightcurves follow as 



F 



F 



(a) 



Cf ' 



F 



00 



F 



00 



F 



Cf 



while the baseline fluxes F,' s ' and the blend ratios o^' are deter- 

basc 

mined as 



(83) 
(84) 
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00 
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Cf 
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F 



00 
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(85) 



(86) 



so that Fg , F B , F^^ e and g^- 1 depend on no other parame- 
ters than ("f, Af, and the caustic passage flux parameters F-' s ' and 
Ff^ s \ where Fg S ' ) even depends on and F-' s ' only. 

The choice of the caustic crossing angle cj> determines the time 
t+ = £e P* in which the source moves by its angular radius (9* 
together with the fold-caustic model parameter as 



tf 



(87) 



For given (d, q, £) determining y { and rif, <f> fixes the location of 
the source trajectory determining the parameter un. 



(VA) f , 




caustic 



Figure 4. Definition of the angles 0, Tf , and r t which describe the relative 
orientation between the source trajectory characterized by ± y f , the caus- 
tic characterized by its inside-pointing normal nf, and the gradient of the 
magnification (VA)f of the non-critical images at the fold singularity y f . 
Regions are marked for which the source trajectory matches the indicated 
conditions: sin <f> > (0 < <f> < 7t) follows from the requirement that 
± j/ f points inside, while the condition Q > requires either cos Tt > 
for u)t > 0, cos r t < for < 0, or cos T t = for uil = 0. 



The time-scale of motion ts = t& sin cf> ^ f B depends both 
on Q and (j>, and with Ff reads 



ts = -75- sin 
Hf 



(88) 



Finally, to is determined by tf = tf ± t * , the location of the 
source trajectory given by y f , nf and <j>, and £e depending on ff, </> 
and Ff . 

Table □ summarizes the dependencies of the quantities char- 
acterizing the microlensing event, namely the proper motion fj, and 
the time-scale of motion 4b as well as their transverse compo- 
nents fi ± and f E , the time i* in which the sources moves by its 
angular radius 9* relative to the lens, the source size parameter 

= /#e = and the source and background fluxes Fg 

and Fg S ', on the fold-caustic parameters (t+, F r , F f * , tDf ), the 
parameters of the singularity (Rf, nt, Af, ( V^4)f ), and the adopted 
caustic crossing angle <f) and caustic rise parameter ft . 

7.3 Parameter search using ujf 

If the complete set of 3 + 2n fold-caustic model parameters 



(tf,t+, Ft s ' , Ff w , uif ) is used in the search of the remaining pa- 
rameters required to characterize the full lightcurve, the parameter 
uif provides a relation between the caustic crossing angle <f) and the 
caustic rise time ft, so that the parameter search is reduced to the 
four-dimensional subspaces (d, q, £, <j>) or (d, q, £, Q). 

As an alternative to expanding ^4Jj thcr (t) which is assumed to 
vary slowly during the progress of the caustic passage around tf, it 
may also be approximated by a linear expansion around tf which 
reads 

A p othet (t)^A f + (t-tf)A { , (89) 
where At = ^othor^f) ( as defined in Sect. |4j and At — 
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Table 1. Determination of quantities characterizing the microlensing event from the fold-caustic parameters (t+ 
singularity (Rf, rif, Af, (VA)f), and the adopted caustic crossing angle <f> and caustic rise parameter Cf. 



F} s \ F^ s \ &t), the parameters of the 



parameter 


caustic 
passage 


singularity 


crossing 


caustic 


other 


f-'* 


** 








0* 


tjr sin^ 


t* 










sin <p 


** 










p — Rf t 1 - 

tx _ £i 

C 2 i 

*E = sin0 
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The choice of a binary lens (d, q) and of the singularity at y { characterized by I yields Rf, rif , Af, and ( VA)f . If one disregards Cit, At = A{ and the value 
of ( VA)f is redundant, whereas otherwise A* = Af — Cf <S? depends on Cf and uit, and u>*, Rf, rif and (VA)f constrain <j> and provide a relation to C,f. 



^other(*f)' while the corresponding expansion of A^ thel (y(t)) in 
space around y f reads (cf. Eq. <44t ) 

^5th«(v(*)) - ^ + (* - **) Vf ■ (VA) f . (90) 
The comparison of Eq. 1891 with the expansion around if, Eq. I64i . 
yields the relations Af — Af and Af = Af =p Af , and by com- 
paring Eq. 490J with Eq. J89}, one obtains 

At = y f ■ (VA) f 

= ±|y f ||(VA) f | cosr t , (91) 

where T t denotes the angle measured from ±y f to (VA)f. As a 
function of the angle ^ Tf < 2n measured from the caustic 
normal rif to (VA)f, which is determined by the properties of the 
lens mapping at y t alone, and the caustic crossing angle < (f> < 
7t, r t reads 

r f 



r t 



-271 



27t 



+ 



(92) 



forcing - f < T t < ^ . With this definition, < cf) < n in- 
creases with a counterclockwise rotation of ± y t . The definition of 
the angles <j>, Tf, and T t is illustrated in Fig.0 

With ("f > and Af = ±Cf ^f > Eq. <9H requires the sign 
of cosTt and uif to coincide. With Eq. 1921 and the condition 
< < 7t, this results in the allowed ranges for Tt and <j> shown 
in Table|2| Configurations (d, q, £) for which there is no viable tra- 
jectory have to be discarded. 

By inserting 



cos Tt = cos Tf sin <fi — sin Tf cos < 



(93) 



and 

(if* : 
LOf : 



1 2/f | = 1/tE as given by Eq. <88> into Eq. <9H . and using 
= ±Af /Cf, one obtains 



flr |(VA) f 



Kt{<t>), 



(94) 



C f 3 * 
where 

iff ((/>) = cos T f - sin T f cot <j) . (95) 

For <2>f 7^ and sin Tf 7^ 0, this yields a relation between 4> and Cf , 
and, according to the previous subsection also between the time- 



scale t-E and the source and background fluxes F s 
spectively. With 



(<0 



and fS. 3 ' , re 



ft|(VA) f | 



1/3 



one finds explicitly 

Cf(0) = Cf,rcf [Kf^)] 1 ^ 



Mf\ 71 , + (Cf/Cf,ref) 3 ~ 

<P«!) = o + arctan Zr-rT- 

Z sm I f 



cos r f 



(96) 



(97) 



(98) 



For Ff = 0, adopts the value Cf.ref (while Uf > 0), whereas 
— Cf.ref is adopted for Tf = 7t (while ujf < 0), irrespective of 4>. 
For uif = 0, one finds that <f> = Tf for T f G (0, 7t) or </> = Tf — 7t 
for Tf € (7t, 27t) . Therefore, the remaining parameter space can 
be parametrized by (d, q, £, <f>) for u>f 7^ and by (d, q, £, Cf ) for 
sin Tf 7^ 0, whereas otherwise such a parametrization is not viable 
due to the adoption of a single fixed value. 

In analogy to Cf, one obtains Af — ±Cf£>f as function of the 
caustic crossing angle <f> as 



AK4>) 

where 

Af, TC f - 



1/3 



(d,n 2 i? f |(vA)fi 



1/3 



(99) 



(100) 



The function [KtiW = Cf (<£)/Cf,ref = A* f ((f)/ A^ tei is 
shown in Fig.[5]for selected angles Ft . [Kf ((j))] 1 is strictly mono- 
tonic in 4> except for Tf = and Tf = 7t, where the constant value 
±1 is adopted. With the requirement of Cf (0) > and the sign of 
Af (cj>) to match that of uif, the range of (j> is restricted as shown 
in Table H For <p -* or <j) -> 7t, | [Kf (4>)] 1/:i tends to infinity, 
whereas zero is approached for (f> — » Tf or (f> — > Tf — 7t, respec- 
tively. 

With the caustic rise parameter Cf depending on the caustic 
crossing angle <f> as described by Eq. <97t . the time-scale of trans- 



verse motion t E and the source size parameter p* 
become functions of d> which read 



m<t>) 

and 



t E , Icl [K ( (<P)] 2/3 



ti/tv also 



(101) 
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Table 2. Constraints on the orientation of the source trajectory imposed by the sign of ui* depending on the orientation of ( VA)f with respect to the caustic. 
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no viable trajectory 



The fold-caustic model parameter lj* characterizes the temporal variation of the magnification due to non-critical images, Tf = /.(rif, (VA) f ) < 2tt 
denotes the orientation of its gradient (VA)f relative to the caustic inside normal rif at the fold singularity y t , I\ = Z(±j/f , (VA)f) is the angle between 
the inside pointing tangent to the source trajectory ±y f and (VA)f, and <f> denotes the caustic crossing angle. 



s 

-e- 




Figure 5. Variation of the magnification of the non-critical images A* as 
well as the caustic rise parameter Q as a function of the caustic crossing 
angle <j> for selected values of the angle Tf between caustic normal rif and 
gradient of the magnification ( VA)f of the non-critical images at the fold 
singularity y f . With . given by Eq. 1 1001 and Cf ,rcf given by Eq. 1961 . 
the curves show A* W/i* ref = Cf/Cf.ref = [#fW>)] 1/3 for T £ = 
0, tt/32, tt/16, ti/8, ti/4, tt/2, 3tt/4, tt, and 3tt/2. The range for the 
crossing angle <\> is restricted as shown in Table l2l so that Q (<f>) > and 
the sign of A* (0) matches that of Cj* . 



P*(0) = P*,re( [A-f(0)n 2/3 

respectively, where 

t B , f e f =W l(VA)f| 

and 



2/:.i 



1/3 



2/3 ti 



(102) 



(103) 



(104) 



For r f = or T f = tt (i.e. sinTf = 0), t£ = rof 



and 



tt/32 





n/8 


k = 






U/4 / \ / 






\ / \tt/2 / 






\ / \3tt/2 / 



'37T/4 



3tt/4 



tt/32 



_L 



tt/4 



3tt/4 



p*, rc f for any 0, whereas otherwise t E (0) and p-i(<f>) are 



Figure 6. The dependence of the time-scale of transverse motion f ^ and 
the source size parameter p* = (5* /(Jg = /tj^ on the caustic crossing 
angle <j>. The upper panel shows curves for ((f>)/t^ t re f = [iff (</>)] 2//3 
while the lower panel shows the corresponding curves for p* (</>) / p* re f = 
[iff (</>)] ~ 2 / 3 for the angle between caustic normal rif and gradient of the 
magnification of the non-critical images at the fold singularity j/ f assuming 
the values T f = 0, tt/32, tt/16, tt/8, tt/4, tt/2, 3tt/4, tt, and 3tt/2. 
The reference values t-g, rc f an <i P* ref me given by Eqs. fKBl and IToIl . 
respectively. According to the sign of il>^ , the caustic crossing angle tf> is 
either restricted to the range Tf < <f> < 7tor < <fi < Tf for Tf G (0, tt) 
or to Tf — tt < < TTor < < Tf — TTf'or Tf G (tt, 2tt) as shown in 
Table|3] so that [K f (<f>)] 2 / 3 and [iff (cj>)] _2 / 3 are strictly monotonic within 
the allowed range except for Tf = or Tf = tt, where a constant value is 
adopted. 



strictly monotonic within the allowed parameter range for <f> shown 
in Table [2] where reaches zero and tends to infinity for 
(f> — > Tf or (f> — *• Tf — 7t, while tends to infinity and p* ap- 
proaches zero for 4> —> or (f> — *• 7t. 

The time-scale of motion £e can also be written as function of 
the caustic crossing angle <f> reading 

t E (0)=t E ,ref [Kt(<t>)] 2ri sin0, (105) 
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Figure 7. The time-scale of motion t E as a function of the caustic crossing 
angle <j> for different selected angles Tf between caustic normal rif and 
gradient of the magnification (VA)f of the non-critical images at the fold 
singularity y t . The curves show t E (</>)Ae rof = [-^"f(0)] 2 ^ 3 sin</>,where 
*E,rcf is given by Eq. 11031 . for T { = 0, n/32, tt/16, tt/8, tt/4, 7t/2, 
37i/4, 7t, and 37t/2. According to Table|2] the condition (f > restricts 
(j> either to the range Tf < <f> < n or < (j> < T[ for Tf e (0, 7t) or to 
rf — 7t< </> < 7tor < </> < Tf — 7tfor Tf 6 (7t, 2tc) depending on the 
sign of iDj . Therefore, t E (</>) nas a unique maximum within the allowed 
range of caustic crossing angles (j>, whereas ts(<j>) approaches zero at its 
boundaries. 



with tE,rcf as defined by Eq. I103i . so that i E reaches a unique 
maximum within the allowed range for (f> and approaches zero at 
its boundaries, located at <j> — 0, (f> = 7t, (f> = Tf or <j> — Ff — it. 
Antiparallel orientations of the gradient of the magnification of the 
non-critical images relative to the caustic normal for which Ff dif- 
fers by (an odd multiple of) n yield identical values for i E , p*, and 
t E . Fig.0shows£^(0)/i E ,rcf = [Kt{</>)] 2/3 andp*(0)/p*, rcf = 
[ivfO)]- 2/3 while Fig.0showste(<?i)/t E , rcf = [K { (<f>)} 2/3 \ sin^ 
for the same selected values of Ff as chosen for Fig. [5] 

Depending on the sign of u)f , A^ tbei (t) as given by the expan- 
sion in Eq. 1641 yields unphysical values A^ theT (t) ^ 1 for either 
sufficiently large or small t. For the region containing data near the 
caustic passage and at least during the caustic passage itself, one 
should require A^ thcl (t) > 1. This condition yields an upper limit 
for the caustic rise parameter (f , which is limited to the range 

Cf < Cf.var = -prjj— (106) 

where r ^ ti denotes the timespan from tf toward the outside 
of the caustic for tl> f * > or toward the inside of the caustic for 
d)f < 0, so that the constraint arises from fulfulling ^ thcr (t) > 1 
at the time tf — t for ± uif > and U + r for ± Ci)f < 0, while there 
is no constraint for tUf — 0. By means of Eq. I98t . the constraint 
on ("f replaces the boundary at <j) = or cj> — tt with </) V ar = <KCf)> 
where the allowed range depends on the sign of &f and the orienta- 
tion of (V4)f with respect to the caustic as shown in Table|5| For 
Tf = or Tf = 7t, Cf adopts the fixed value |Cf,ref| irrespective 
of the caustic crossing angle <f>, so that no restriction on <j> results if 
ICf.ref | < Cf,var, whereas otherwise the chosen configuration char- 
acterized by (d, q, £, Cf ) has to be discarded. 



7.4 Blending constraints 



In the absence of blending, i.e. = or g {s) — for all 

lightcurves, Eq. 1701 fixes Cf to 

A[ 



(107) 



where g f * is a fold-caustic model parameter. If one disregards (2> f * , 
the parameter search is reduced to the four-dimensional subspace 
(d, q, £, <j)). Otherwise, Af — Af — Cf u)f ti, so that 



Cf = 



Af 



g\ + ti 



(108) 



and <f) is determined by means of Eq. I98i . leaving only the three- 
dimensional subspace (d, q, I) to be explored. 

Restricting the parameter space to non-negative background 
fluxes, i.e. ^ or ^ 0, yields a lower limit to Q, limiting 
it to the range 



Cf > Cf,B = 



At 



(109) 



with gf min > being defined by Eq. J73L The condition gf min > 
corresponds to the fluxes for all sites at time tf (when the leading 
limb of the source enters or its trailing limb exits the caustic) to 
fulfill F* (a) > for F^ s> > 0, whereas at least one of the con- 
ditions F^ s) > and Af > 1 would be violated if F* (s) < 0. 
Similarly, the simulataneous fulfillment of F^ > 0, Af > 1 and 
Fjj implies the fluxes at time tf (when the source center crosses 
the caustic) to be restricted to Ffoid(if) = F*' s ' + Ft" uif t+ 
= F^ s) A{ + F^ } > 0, so that ^ min + d) f * ti > needs to be 
fulfilled, posing a further restriction in the case cD f * < 0. 

The blending limit on Q given by Eq. J109I implies a lower 
limit to the time-scale of perpendicular motion 



r 2 

■ X > t.f,B f 

tE> lif- t > 



and an upper limit to the source size parameter 



Rt ti 



(110) 



(HI) 



while the ranges of the fluxes and blend ratios are restricted to 








r s 


^B 
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^ r (s) /Cf 



F^gt 
p(») 



9f : , 



0. 



By definition, 



F, 



9i, 



F 



(112) 
(113) 

(114) 
(115) 

(116) 



where the equality holds for the lightcurve s mln involving the min- 
imal value g*f ( " min) = 3 f * min . 

With 0(Cf) given by Eq. 1981 if one takes the fold-caustic 
model parameter uif into account, the lower limit on due to non- 
negative backgrounds translates into a limit </>(Cf ,b) = 0(Cf,s) for 
the caustic crossing angle (f>, which replaces the boundary at Ff (for 
T f e (0,7t))orr f -7t(forr f e (n,2n)) as shown in Table[3] With 
Cf = ICf.rcf | for rf = or F{ = n irrespective of the caustic cross- 
ing angle (j>, no restriction on <f> results if |Cf,ref| ^ Cf.B, whereas 
the configuration (d, q, £, provides no viable solution otherwise. 
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Table 3. Constraints on the caustic crossing angle <j> due to a limit on the variation of the magnification of non-critical images and the requirement of non- 
negative background fluxes. 
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The upper limit Q < Cf.var which avoids A p , (t) < 1 over the temporal range covered by observations and the lower limit Q Cf ,B which avoids 
negative background fluxes translate to the shown limits on the caustic crossing angle </> V ar = <fi(Q ) or <Pb = '/'(Cf b)> respectively, with <f)(Cl ) given by 
Eq. 1981 . where the allowed range depends on the temporal variation of the magnification due to the non-critical images described by the fold-caustic model 
parameter Cjf and the the orientation of their gradient (V.A)f given by the angle Tf between the caustic inside normal n.f and (V-A)f. The variation limit 
4>va,T replaces the boundary at <fi = or tf> = Tt, whereas the blending limit 4>B replaces the boundary at <fi = Tf or (f> = Yf — Tt. 



7.5 Constraints from data outside the caustic-passage region 

The tails of the lightcurve of a binary lens microlensing event ap- 
proach those of a single lens event. From data taken in the cor- 
responding regions (far outside the caustic passages), the baseline 
fluxes F^jL, are easily determined while data of better quality will 
allow to determine the time-scale of motion tE and the blend ratios 



g (s> , yielding measurements of the source fluxes Fg S ' and the back- 
ground fluxes F B . As will be pointed out later, measurements of 
these parameters are vital for constraining the mass ratio q and the 
angular separation d 8e of the binary lens and for predicting other 
caustic passages. 

If the source flux Fg is known for (at least) one site k, the 
caustic rise parameter (f is determined as 



Cf 



F 



Ik) 



(117) 



while it follows from the background flux F^\ baseline flux F^ e 



or blend ratio g^ k ' with Af as 



Cf = 



Cf 



F, 



(fc) 



F, 



(k) 



F r W Cj* t± - F B k) 



F, 



(fc) 



rp*(k) , rp(k) j_ rp(k) 



Cf 



F 



(k) 



f; 



(k) 



(A { + 



{A f - 1) 



(118) 



(119) 



(120) 



F W _ F > 



The conditions F s (fe) > and At > 1 imply F fo u(t{) = F s (fc) A f + 
(fe) + F r (fc) Colli > F B k) and F* (k) + F r (fc) Qfti > 
F blL> so that Cf > is ensured. If F* (k) + F r (fe) ii f * = 0, Cf 
cannot be determined from g^ k \ whereas Af follows directly as 
A{ — — <r ' in this case. 

The coefficient Af = Af — Q cDf t+ in general depends on the 
fold-caustic model parameters Cjf and and on the choice of Cf, 
while for ti f * = 0, this relation reduces to Af = A{ and Af follows 
directly from the choice of (d, q, £). 



fluxes and the relative blending between the different lightcurves, 
while their absolute values are related to (f and Af. Explicitly, 
Eq. 1651 implies that the source flux for any lightcurve s is deter- 



mined by the source flux for a specific lightcurve k with F r and 

F, 



(fc) 



F W - zi F 



F 



(fc) 



(121) 



while Eqs. IS 1 1 and 1821 yield the background and source fluxes 
for any lightcurve from the corresponding values of a specific 
lightcurve with F r (s) , F r (fc) , F* (s) , 



and F, 



:(fc) 



F 



?(.*) 



F 



= F, 



:(,) 



F 
F 

F (s) 
F, 



( fe) ^f 
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Fr S) ( 



(fc) 



(fc) 



F 



(fc) 



(122) 



(123) 



By means of Eqs. 11 171 to 1 1 201 . Af can be determined if two 
of the parameters Fg k \ F B k \ F^* s ' c , or </ fc ' are known for the same 
lightcurve, which is already the case if two of these parameters are 
known for any pair of different lightcurves. 

Similar to the fluxes, the blend ratios $' s ' for all sites can be 
deduced from the value for any site g^ k ' as 



9 



Af 



Af 



(124) 



but in contrast, this relation involves the parameter Af, so that the 
measurement of two different blend ratios </ fe ' ^ ^(0 yields not 
only Q as a function of Af, but also 



Af = 



F} k) {f; 
Fj l) (f; 



(fc) + F (fc) ». 



Wf t* ) 



(125) 



With the values of the caustic rise fluxes F r , the determina- .. I(T , .... (k) , (n . t . . . 

n ; l r t ^ i itself The condition (p ' ^ g l ' ensures that the denominator does 



tion of Cf fixes all source fluxes Fg according to Eq. <831 . and to- 
gether with the caustic background fluxes F f * ^ and Af, all back- 
ground fluxes Fg S ', baseline fluxes F^^ e and the blend ratios 
follow according to Eqs. 1841 to 1861 . Therefore, a measurement 
of one of these quantities for a single lightcurve fixes all of these 
quantities for the full set of lightcurves. 

Is) * (s) 

The fluxes F r and F f obtained from modelling the data 
near a caustic passage contain information about the relative source 



not vanish. For Fl 



+ F r w w? ti 



-g 



(fc) 



0, Eq. 11251 reveals At 



making the measurement of g"' redundant for obtaining Af 
in this case. 

For a given (d, q,£), the determination of Cf not only yields 
all fluxes and blend ratios, but also the source size parameter 
and the time-scale of perpendicular motion i B . If uif is disre- 
garded, the parameter search reduces to the four-dimensional sub- 
space (d, q, £, <j>). Otherwise, with a measured ujf, Qt determines cj> 



Theory and practice of microlensing lightcurves around fold singularities 15 



with Rf , nf and (VA){, so that the proper motion /i, the time i* in 
which the source moves by its angular radius, the time-scale of mo- 
tion £e, and the whole lightcurve are defined, while the search of 
parameters is reduced to the three-dimensional subspace (d, q, £). 

A configuration (d, q, £) can be discarded if Ai does not match 
the value determined from a measurement of two fluxes or two dif- 
ferent blend ratios or from a flux and a blend ratio. 

A measurement of the time-scale of motion £e yields a relation 
between the caustic crossing angle <f) and the caustic rise parameter 
Cf , depending on Rf defined by (d, q, £) and the fold-caustic model 
parameter tf, thereby reducing the number of free parameters by 
one. 

The position of another caustic passage provides a relation be- 
tween t-E and (f> which also implies a relation between Cf and <f) with 
Ri. 

Since each of the measurements of Uf , tE, and the position 
of another caustic passage provide a relation between Cf and cj> de- 
pending on Rf, rif and (Vi)f, and a measurement of one of the 
fluxes Fc! s) , F^ s) , F b ( ^ e or a blend ratio g (k) yields Cf, both < f and 
(j> are determined and the space of free parameters reduces to the 
three-dimensional space (d, q, £) with two of these parameters be- 
ing determined, while the determination of more than two of these 
parameters yields constraints for (d, q, t) which characterize the bi- 
nary lens model and the position of the singularity on the caustic. 
With the restriction of non-negative background fluxes Fg S ' ^ 
and Cf being determined, Eq. i 1091 yields a limit 

At < Cf (at. + «f **) i ( 126 > 

which puts a restriction on the choice of the binary lens and fold 
singularity parameters (d, q, £) for a given caustic crossing angle 

<t>. 

The arising constraints for (d, q, £) and the power of the pre- 
diction of other caustic passages are further discussed in the fol- 
lowing section. 



8 PREDICTIVE POWER 

A dense photometric sampling of binary lens microlensing events 
that involve fold-caustic passages provides opportunities for resolv- 
ing the stellar atmosphere of the source star yielding measurements 
of limb-darkening coefficients, for measuring the proper motion fi 
of the source relative to the lens, and for obtaining the physical 
properties of the lens system such as the total mass, mass ratio, 
semimajor axis, and orbital period. From all lens properties, only 
the mass ratio q and the separation parameter d — 8/9e are observ- 
able in the lightcurve, whereas the lens mass M and the distance 
Dl are convolved into the time-scale of motion t-E which however 
can be disentangled from an assessment of effects in the lightcurve 
caused by the parallactic motion. The measurement of semimajor 
axis and orbital period is hampered by the fact that the lightcurve 
depends on the instantaneous angular separation only (except for 
some small effects caused by orbital motion) leaving inclination 
and eccentricity of the binary system unconstrained, so that these 
quantities can only be assessed statistically. 

The determination of the local properties of the lens mapping 
near the fold singularity parametrized as described in Sect.|6|allows 
to study the atmosphere of the observed source star and to obtain a 
lower limit on its proper motion /i relative to the caustic (i.e. rela- 
tive to the lens system approximately), namely its component fx 
perpendicular to the caustic. 

A powerful test of stellar atmosphere models requires a 



dense coverage of th e lightcurve during the caustic passage (e.g. 
iRhie&Benneti 19991) . In order to be able to schedule the necessary 
observations, including taking some high-resolution spectra which 
will provide a deep probe of the chemical composition of the source 
star by means of observable variations in associated spectral lines, 
caustic passages need to be predicted some time ahead. Unless the 
source trajectory approaches a cusp, fold-caustic passages appear 
in pairs, comprised of a caustic entry and a subsequent caustic exit. 
Caustic entries are practically unpredictable before they actually 
occur, so that their beginning phase can only be caught by chance 
and immediate action has to be taken to follow the lightcurve over 
the caustic entry. Since the behaviour of the lightcurve around a 
caustic passage depends on local properties only, the data taken dur- 
ing the caustic entry contains practically no information about the 
follow ing caustic exit (c.f. lAlbrow et alll999btlJaroszvnskT& Mao 
l200ll) . However, caustic exits can be predicted using the data on the 
rise to the caustic passage peak. Once such a rise has progressed, 
the parameters t| = tf, Fj s \ and Ff can be determined for a 
point source (t+ = 0) and = 0. When the source reaches the 
region where the curvature of the lightcurve changes sign, it be- 
comes possible to assess the source size by including t* > in the 
fit, thereby obtaining a prediction for the end of the caustic exit at 
tf which however depends on the amount of limb darkening which 
might be determined from a fit to the data as well as u)f at later 
stages. 

The data taken near the caustic passage also neither provide 
constraints on the mass ratio q and the angular separation parame- 
ter d = S/8e which characterize the binary lens nor on the location 
of the fold singularity y f on the corresponding caustic described 
by the parameter I. The parameters (d, q, £) are related to the fold- 
caustic model parameters by means of the local properties Rf, nt, 
Af and (VA)f which are not observables themselves but affect the 
lightcurve through convolutions with the source size parameter 
and the source and background fluxes and F^K The insen- 
sitivity to source motion parallel to the caustic leaves the caustic 
crossing angle <j> undetermined. The parameters (d, q, £, <fi) fix the 
binary lens system, the spatial position of the source trajectory, and 
with tf — tf ±f* the point of time where the source center crosses. 
The unknown source flux (or baseline flux) results in a free caustic 
rise parameter Cf- With the rise fluxes _F r , Cf yields the source 
fluxes Fg a \ and together with the caustic baseline fluxes Ff^ 
and Af, 12 it yields the background fluxes F^ . Together with Rf, 
Cf yields the source size parameter = O^/Oe and the time-scale 
of transverse motion tg, while the time-scale of motion t-E also 
depends on <f>. 

Since the observed caustic-passage time is the product of 
the transverse time-scale of motion tjjr =tn/ (sin (f>) and the source 
size parameter p+, smaller sources with smaller transverse proper 
motion or larger sources with larger transverse proper motion can 
yield the same lightcurve. Moreover, the insensitivity to the cross- 
ing angle means that the same transverse proper motion can be 
achieved for any angle by a corresponding choice of the full proper 
motion. 

While the determination of the model parameters for the 
lightcurve near the caustic passage alone cannot yield any predic- 
tions for its remaining parts, the combination with a few simple 
characteristics of the data outside the caustic-passage region how- 
ever provides some vital constraints. 



Aj = Af — Cf uif t-^r , which reduces to = Af for Clf = 0. 
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Restrictions of the parameters (d, q, £) can arise from indepen- 
dent measurements of (combinations of) the local properties Rf, 
fit, At or (VA)f. As pointed out in Sect, 17.51 measurements of 
one baseline flux and the blend ratio g' fe ' or of two blend ra- 

basc •J 

tios g' ^ and g"' yield Af, iz while measurements of Cjf, one base- 
line flux iv?2 e , an d the time-scale of motion provide a relation 
between Rf, At nt, and (VA)[. 

For fixed fold-caustic parameters (tjr, Fj s \ , ilif), the 

full lightcurve is defined with the choice of (d, q, I, Q , (f>), so that 
all of its features such as the position, duration and peak flux of 
other caustic passages or other types of peaks, and the time-scale 
of motion <e are determined. 

By requiring a match of the model lightcurve with the ob- 
served data, the determination of Q and <j> or the provision of a rela- 
tion between these parameters together with limits, arising from the 
appropriate sign of Af (that allows to fulfill ff > with the sign of 
u)f) or the restriction to positive background fluxes, therefore yields 
restrictions for the choice of (d, q, I). Alternatively, the restrictions 
on (d, q, I, £f , (f>) constrain the behaviour of the lightcurve outside 
the caustic-passage region, so that in particular limits on the posi- 
tion, duration, and peak flux of other caustic passages arise. Due to 
intrinsic ambiguities for binary lens models iDominik 1999 albl) . a 
unique solution however cannot be expected and different scenarios 
will be possible. 

A relation between ff and (j> is provided by a measurement of 
(if (which also yields an upper limit on £e) or of tv, and together 
with a measurement of a baseline flux F^ e > Cf an d <j> are deter- 
mined for given (d, q, I). 

Determining Cjf and t% fixes C,t and <j) with Rf, rif, and 
(V.A)f. If a baseline flux F^l B is known in addition, a relation 
between R{, nt, At, and (VA)f results. 

The observed position of another caustic passage provides a 
relation between £e and (j> which implies a relation between Q and 
4> with Rf. Combining this information with measurements of one 
of the parameters Qf , ts or a baseline flux F^ e f* xes Cf an d <^> 
while a combination with more than one of these parameters also 
yields relations between Rf, rif, Af, and (VA){, thereby restrict- 
ing (d, q,£). Other features of the additional caustic passage such 
as its duration, its peak flux, or the flux and rate of change of flux 
at the limb of the source entering or exiting the caustic, yield addi- 
tional restrictions which can become quite powerful (cf. lAn et alJ 
l2002h . 

If £f is determined, the requirement of non-negative back- 
ground flux Fg S ' ^ constrains Af, so that the choice of (d, q, £) 
is restricted if <j> is determined in addition. 

For a worked example, some considerations about the power 
of different characteristics of the photometric data outside the 
caustic-passage region for determining the lens properties and pre- 
dicting other caustic passages have been made by Albr ow et alJ 
l!999bh . 



9 SUMMARY 

Dense high-precision photometric sampling of microlensing 
lightcurves during fold-caustic passages can provide both a mea- 

13 -'base ' s easu y obtained as observable of the lightcurve, while is 

(k) 

also determined with F, and a measurement of either the source flux 

base 

or the background flux F^\ 



surement of the proper motion p between lens and source star and a 
probe of the stellar atmosphere of the source star, e.g. parametrized 
by limb-darkening coefficients, allowing to test existing theoretical 
models. Making use of the alerts supplied by the current microlens- 
ing surveys OGLE and MOA, a follow-up campaign like PLANET 
is capable of providing 10-15 such measurements per year on dif- 
ferent types of galactic stars. 

While the duration of the caustic passage tj" obtained from 
modelling the data in the caustic-passage region and of the angu- 
lar size of the source which can be obtained by spectral typing, 
directly yield the proper motion p x = of the source perpen- 

dicular to the caustic, the full proper motion p = p x / (sin (f>) only 
follows with a determination of the caustic crossing angle (j> from a 
model of the full lightcurve. 

During fold-caustic passages, the microlensing lightcurve 
shows an increased sensitivity to the parallactic motion of the 
Earth around the Sun and to the orbital motion of the binary lens 
whose measurement breaks the degeneracy that remains between 
lens mass M and distance Dl a fter determination of the proper 
motion p iDominik 1998; An et al. 2002). Therefore, microlensing 
events that involve caustic passages are prime objects for determin- 
ing the mass function and the phase-space distribution of stellar 
populations whose constituents act as gravitational microlenses. 

Moreover, the prominent feature of caustic-passage peaks al- 
lows a proper determination of binary lens parameters and with the 
possibility to measure p, M and Dl individually, a reasonable op- 
portunity is provided for obtaining statistical distributions of prop- 
erties of stellar and sub-stellar binaries such as their mass, mass 
ratio, semimajor axis, and orbital period. 

For a given distribution of lens masses, the lens mapping maps 
the position of an observed image of a source star to its true posi- 
tion. While there is a unique true source position for any observed 
image, a source may have several images. While the lens mapping 
is regular (i.e. its Jacobian determinant does not vanish), it can be 
locally inverted, so that the number of images only changes if the 
source reaches a singular point of the lens mapping. The lowest- 
order singularities are folds which form closed curves called criti- 
cal curves in the space of image positions and caustics in the space 
of true source positions. When a point source crosses a fold caustic 
from the 'inside' to the 'outside', the magnification of two of its 
images diverges with the inverse square-root of the distance per- 
pendicular to the caustic, where these two images become criti- 
cal, while they disappear just as the source exits the caustic. The 
fold caustic is characterized by the single parameter R{ describing 
its strength, i.e. a proportionality factor between the magnification 
and the inverse square-root of the perpendicular source-caustic dis- 
tance. 

The magnification of the critical images of an extended source 
with source size parameter p* = 9±/9e is characterized by a uni- 
versal caustic profile function describing the shape of the lightcurve 
during a caustic passage which is characteristic for any given radial 
brightness profile of the source, while for the same caustic strength 

— 1/2 

R{, the magnification is proportional to p* . If the luminosity of 
the source at the limb does not vanish, a discontinuous change in 
the slope of the lightcurve is produced as the leading limb of the 
caustic hits the caustic. On the source further entering the caustic, 
the magnification rises to a peak and then drops asympotically pro- 
portional to the inverse square-root of the perpendicular distance of 
the source center from the caustic. 

The caustic-passage regions of n observed lightcurves (for dif- 
ferent sites and/or spectral bands) during a caustic passage can be 
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characterized by the 3 + 2n parameters (tf , t^, F T • , Ft , U)f). 
The source crosses the caustic during the timespan 2 t t , where its 
leading limb enters or its trailing limb exits the caustic at time tf 
where the flux is given by F { * . Fx denotes the flux of the two 
critical images at a unit time t after the source center enters or be- 
fore it exits the caustic that are produced if the extended source is 
replaced by a point source at its center, and the parameter Qf mea- 
sures the rate of the (constant) change of flux in units of iy due 
to the non-critical images for the source moving inwards, while this 
rate is given by — u)j! for a source moving outwards. 

Since the influence of the lens mapping on the microlensing 
lightcurve near a fold-caustic can be approximated by means of 
local properties, the study of stellar atmospheres, and in particu- 
lar the determination of limb-darkening coefficients, does not re- 
quire the assessment of a complete set of model parameters that 
characterize the full light curve and is not influenced by ambi- 
guiti es and d egeneracies th at are likely to occur in these pa rame- 
ters jp ominik & Hirshfeld 1996; lDominikiri999albt lAlbrow et all 
1999b). However, this also means that the data taken in the caustic- 
passage region does not provide information about the behaviour of 
the lightcurve outside this region. In particular, the caustic-passage 
data do not allow to obtain predictions about subsequent other caus- 
tic passages. 

For a binary lens model and a fold singularity on its caus- 
tic, all local properties of the lens mapping can be obtained an- 
alytically by means of derivatives of the Fermat potential. How- 
ever, these properties are not observables themselves but affect the 
lightcurve through convolutions with the source size parameter p* 

( s) Is) 

and the source and background fluxes Fg and F^ which in gen- 
eral are unkown themselves. The data in the caustic-passage region 
therefore neither constrains the mass ratio q of the binary lens and 
the angular separation S = d #e between its components, nor the 
location of the fold singularity on the corresponding caustic. 

Besides parameters that characterize the source brightness 
profile, a microlensing event involving an extended source and a bi- 
nary lens is characterized by 7 + 2n parameters. By modelling the 
data around an observed fold-caustic passage, thereby fixing 3 + 2n 
parameters, or 2 + 2n parameters if Ci)f is disregarded, the search 
for model parameters describing the full lightcurve reduces to a 
search in a four- or five-dimensional subspace holding the remain- 
ing parameters which can be chosen as (d, q, £, <j>) or (d, q, £, <j>,Q), 
respectively, where £ characterizes the position of the cusp singu- 
larity on the caustic and ff characterizes the rise of the magnication 
on the approach to the caustic. The free choice of (j> reflects the in- 
sensitivity to source motion parallel to the caustic, while the free 
choice of £f refl ects the freedom in ch oosing the source flux. As 
demonstrated by Albr ow et alj ll999bl) . this approach provides an 
efficient method for obtaining all suitable models that are consistent 
with the data of the observed microlensing event. 

The determination of a few simple characteristics of the data 
outside the caustic-passage regions such as source, background, 
or baseline fluxes, the blend ratio, the time-scale of motion ts, 
or constraints from other characteristic observable features of the 
lightcurve such as the position, duration, or peak flux related to 
other caustic passages or other peaks yields further restrictions on 
the parameter space to be investigated and leads to a reduction of 
its dimension. This information is crucial for obtaining constraints 
on the mass ratio q and the angular separation parameter d = 5/8e 
which characterize the binary lens, and on the properties of subse- 
quent caustic passages. 
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Table 4. Used symbols 



Symbol Meaning 



Definition or prime relation 



Caustic-passage Ughtcurve 

t, 



f 



r 

*(s) 



time of beginning (end) of caustic entry (exit) 
caustic passage half-duration 
temporal variation of non-critical images 
caustic rise flux 
flux at time t£ 
caustic rise flux ratio 
magnitude at time tf 

limb-darkening coefficient for profile oc cos p 1? 
Caustic-passage auxiliary quantities 

tf time when source center passes caustic 

ff caustic rise parameter 

t T caustic rise time 

Binary-lens Ughtcurve 



± f 
9t 

V 



M 

'{ 

M 
M 



*E 

to 
a 
d 
1 
P* 

p s 
F M 

r B 

F (s) 
base 

M 



event time-scale 
lens-source impact parameter 
time of closest lens-source approach 
angle between binary lens axis and source trajectory 
lens separation parameter 
binary lens mass ratio 
source size parameter 
source flux 
background flux 
baseline flux 
blend ratio 
Local lens properties 

R f caustic strength 

TOf caustic inside normal 

y f location of source singularity 

Af magnification due to non-critical images at y { or tf 

Aj: magnification due to non-critical images at tf 

temporal derivative of A^ ther at tj? 
spatial derivative of A^ thcr at y t 
angle from n t to (Vi)f 
Lens-source motion 

<j> caustic crossing angle 

relative lens-source proper motion 
lens-source proper motion perpendicular to caustic 
event time-scale referring to perpendicular motion 
time for the source to move by its angular radius 
y f source proper motion vector per 9 E at tf 

F t angle from ±y f to (VA)f 

Other quantities 

i arbitrary unit time 

M total lens mass 

Mi individual lens masses 

D L lens distance 

Ds source distance 

6+ angular source radius 

S angular separation betwen lens objects 

0n smallest angular lens-source separation 

#E angular Einstein radius 



(VA) f 

r f 



t = P*t E 



r( S ) 



f 



9f ~ -f 

* w 
V 



F^Q 



F^A* 



in,. 



/Ft 

fold 1 - 



t t = t* ±t± 

t r = Rft^ 

u = e /e E 



d = s/e E 

q = M 1 /M 2 



base 

„(*) = 



Vi -■ 

Af 

At 



y(tf) 

z j4 other(ff) 
other 
other 



other > 



.„(**) 

(t*)~ 



(VA) f 
r f = Z(n f ,(VA) f ) 



^ ther (y f ) 



Vl 

r t 



= /XSlIKp 

= t E /(sin<£) 
p*t E = 6*1 'n = ti 
--y(tf),\y t \ = l/t E 
= /(±y f ,(VA) f ) 



M = J2 M * 



9 E = [(4GM/c 2 )(D s - Dl)/(I3s-Dl)] 1/2 



Upper (lower) signs correspond to a caustic entry (exit), F^ 3 X(t) denotes the flux, mfL(t) the magnitude, A p h r (t) the magnification due to non-critical 
images, and y(t) = 9(t)/6 E the source position as a function of time, where 0(t) is the angular position of the source. Among the superscripts and 
subscripts, Of denotes a relation to the fold singularity, □* denotes a reference to the stellar limb, □+ denotes a relation to stellar radius, □ denotes the 
vector component perpendicular to the fold caustic, and d' 3 ' refers to a specific site and filter. 
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